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Abstract The overall properties of the Herbig-Haro 
objects such as centerline velocity , transversal pro- 
file of velocity , flow of mass and energy are explained 
adopting two models for the turbulent jet. The com- 
plex shapes of the Herbig-Haro objects , such as the arc 
in HH34 can be explained introducing the combination 
of different kinematic effects such as velocity behavior 
along the main direction of the jet and the velocity of 
the star in the interstellar medium. The behavior of the 
intensity or brightness of the line of emission is explored 
in three different cases : transversal ID cut , longitu- 
dinal ID cut and 2D map. An analytical explanation 
for the enhancement in intensity or brightness such as 
usually modeled by the bow shock is given by a careful 
analysis of the geometrical properties of the torus. 

Keywords Pre-main sequence objects, young stellar 
objects (YSO's) and protostars (T Tauri stars, Orion 
population, Herbig-Haro objects, Bok globules, bipolar 
outflows, cometary nebulae, etc.; Jets, outflows, and 
bipolar flows 



1 Introduction 

The Herbig-Haro objects , in the following HH , are 
on the short distances from the central source col- 
limated patches of nebulosity associated with a cen- 
tral very young star; the name derives from the first 
two astronomer which studied them in details , see 
Herbigl (|1950|); iHarol (|1952|). The HH are observed 



in the various band of the electromagnetic spectrum 
: in the radio , see Rodriguez and Reipurth (1989); 
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Curiel et al. 


(19891; 1 


\nglada et al 


(19921; Curiel et al. 


(1993 


I;] 


lodriguez and Reipurth' (1994) ; in the infrared 


see Reipurth and Aspin, (1997); , Reipurth, Devine, and 


Bally (19981; 


Chrysostomou et al. 


2000|; 


Davis, Smith, 


and Eisloffel 


(20001; 


Davis et al. 


20021; 


Takami et al. 


(2005 20061; in the optical see Schwartz et al. (19881; 


Rolph, Scarrott, and Wolstencroft (1990); Scarrott 


et al. (1990); Bohm, Raga, and Binette (1991) 


Uchida 


et al. 


( 


1992 


1; Gomez, Whitney, and Wood 


(19981; 


Masciadri and Raga 


(20011; in the ultraviolet 


Dopita, 


Binette, and Schwartz 


(1982); C 


Cameron and Liseau 


(jl990 1; Bohm, Raga, and Binette ( 


1991); 


Boehm, Scott, 


and So' 


f (1 


391); Boehm, Noriega-Crespo, and Solf 


( 1993 1 ; in the X-ray see Pravdo and Angelin 


(19931; 


Raga, Noriega-Crespo, and Velazquez' f2002|). The 



HH's are also observed through emission-line spectra 
. We remember that the astronomers observe both the 
forbidden emission lines from low ionization species , 
[SII] and [OI] as well as ionization lines such as [Oil I] 



and [Nell I], see Hartigan and Morse (2007) 



On long distances from the central source the HH are 



curved, see Salas, Cruz-Gonzalez, and Porras (1998); 



Bally and Reipurth (2001); Bally et al. (2006). From 



a theoretical point of view the apparent deceleration 
of HH34 has been modeled by the interaction of a 
fragmented jet with the surrounding environment, see 



Cabrit and Raga ( 2000 1 . The theoretical problem of the 
curvature of the HH objects has been attached adopt- 
ing an HH jet/counterjet system that is immersed in an 



isotropic stellar wind Raga et al. (2009), or discussing 



the properties and similaxities of the laboratory and as- 



trophysical flows, see Ciardi et al. (20081 and Hartigan 



et al. (2009). The already cited models concerning the 
HHs leave a series of questions unanswered or partially 
answered: 

• Which are the laws of motion that regulate the prop- 
agation of HHs in the Interstellar Medium ? 
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• Is it possible to model the main physical properties 
of HHs such as matter entertainment or mechanical 
luminosity ? 

• Can we model the bending of the HHs at the light of 
the known rules of the kinematics ? 

• Can we model the intensity or brightness behavior 
along the HHs using different models ? 

In order to answer these questions Section [2] reports 
the astronomical data of two HH objects. Section |3] 
reports two physical theories on turbulent jets which 
are converted in astrophysical equations in Section [3j 
Complex phenomena such as jet's bending and train of 
knots are explained in Section |5] adopting the compo- 
sition of different kinematic effects and the theory of 
the Kelvin-Helmholtz instabilities . A set of models 
for the intensity or brightness of HH's which arise from 
geometrical arguments are explored in Section |6] Sec- 
tion [7] analyzes a model for the intensity or brightness 
of HHs as given by a linear and a non linear relationship 
between emission coefficient j and turbulent power. 



2 The astronomical observations 

This section presents the astronomical data of HHl and 
HH34. 

2.1 The inner part of HHl 

The coUimated jet HHl has been observed in different 
astronomical wavelengths such as [Fell] by 



Reipurth 



et at. (20001 , near infrared by Davis et at. 



far infrared by Molinari and Noriega-Crespo 



optical/infrared by Nisini et al. (20051 , Ha by Riera 



Garcia Lopez et al. ( 2008 1 



et al. (20051 , UV by Li et al. (20071 and FeII/H2 by 



20001 



(2002 



HHl is part of the complex HH 1/2 that covers 3' 
with its brightest components. The distance of this 
comple x is 460 pc, see [Molinari and Noriega-Crespo| 
( |2002[ ); [Garcia Lopez et al.\ ( |2008[ ). 

The length of the jet , Lhhi , measured from the 
central source , VLAl , to the knot A is according to 



Nisini et al. (20051 



L 



HHl 



22" = 0.049 D. 



'460 pc , (1) 



ble 2 in Nisini et al. 



see Ta- 

The initial diameter di 
(measured at the knot L-I) and the final diameter df 
(measured at the knot A ) are , according to Table 3 in 



where i346o is the distance in units of 460 pc 
p005| 



Nisini et al. (20051 



= 2 X 0.1" = 0.000446 i:»46o pc 
df = 2x 0.5" = 0.00223 D460 pc ■ 



(2) 



The averaged radius of inner part HHl is r = 0.3" and 
9 blobs characterize the structure. The half opening 
angle , a/2 , is 

^ = arctan '^//^ ~ ^'/^ ^ 0.OI8I rad = 1.041 ° . (3) 

With these data the laboratory parameter a:/(i is is 
X Lhhi 



110 



(4) 



where x represents the jet's length and d the nozzle's 
diameter . 

2.2 HH34 the giant jet 

HH34 constitutes the archetypal bipolar coUimated jet 
from a young star and has been carefully studied in deep 
Ha and [SII] with the Wide Field Planetary Camera 
2 on board of the Hubble Space Telescope and at the 



et al. (20071 



Gemini Observatory, see Reipurth et al. (2002); Beck 



According to the data on HH34 as suggested in Mas 



ciadri et al. (20021 and Reipurth and Raga (19991 the 



distance is 460 pc, the length of the jet (arc comprised) 
, Lhh34 , is 



18 



iffff34 l-5pc or 4.6210 
the initial jet's diameter is 
d^ = 9.7 10~ V or 6 X 10^^ cm 
The ratio x/d, takes the value 



X _ LhH34 

d di 



773 



(5) 



(6) 



(7) 



A detailed study of the inner part of HH34 shows a well 



coUimated jet , see Beck et al. (2007) , which means an 



opening angle of few degree , i.e. a = 2.86° (measure of 



the author). Making reference to Figure 3 in Reipurth 



et al. ( 2002 ) the inner part of HH34 has an averaged 



radius of 0.5", a length of 27" and twelve blobs 



3 The turbulent jet 

This Section reports known and new formulas on tur- 
bulent jets. We selected the modern approach to tur- 



bulent round jets, see Pope (2000) and the classical ap- 



proach of Landau (1987). These two approaches were 
already used to model the extragalactic jets , see [Za-[ 



ninetti (2007 20091. 
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3.1 The exact solution 

The theory of turbulent round jets can be found in dif- 
ferent textbooks. The more important formulas are 



now reviewed as extracted from chapter V in Pope 



(20001 ; similar results can be found in Bird, Stewart 



and LightfootI ( |2002D and in [Schlichting et al.\ ( |2004D . 
We start with the centerline velocity Uo{x) , equa- 
tion (5.6) in iPope (20001 , as measured in the labo- 



ratory experiments : 



Uojx) 



B 

x/d 



(8) 



here x denotes the main direction , d is the diameter 
of the nozzle, _B is a constant derived in the labora- 
tory that takes the value 5.8, and Ui is the initial jet's 
velocity. The solution of the mean velocity < U > , 
equation (5.100) in Pope (20001 , along the main direc- 
tion is 



< 



,d. 1 

U >= Ui8aTiyT{-)j— ^ 

x' (1 + ax'q'^) 



(9) 



where rj — ^ , r is the radius of the jet at x, ax is a 
constant and is the turbulent viscosity. The viscosity 



, equation (5.104) in |Pope| ( [2000[ ) , is 

S 



/2- 1) 



and ar, equation (5.18) in Pope (2000), is 



ax 



(V2~l) 
52 



(10) 



(11) 



where S is connected with the opening angle a through 
the following relationship 



S = tan 



(12) 



The production of turbulent kinetic energy in the 
boundary layer approximation , equation (5.145) in 



Pope (2000) , is 



dy 



(13) 



where y is a Cartesian coordinate that can be identified 
with r, the perpendicular distance from the centerline 
and the units are -^^^^ . The flow rate of mass mix) is 



, see equation (5.68) in Pope (2000) , 



m(x) = 2^p(&i/2(a;))2c/o(x) X / U{i)d^ , (14) 



where 



(15) 



and 

/(0 = 



1 



(16) 



where A is a constant that will be later defined and 
61 is the value of the radius at which the velocity is 
half of the centerline value. The jet draws matter from 
the surrounding mass of fluid. Hence, the mass of fluid 
carried by the jet increases with the distance from the 
source. The flow rate of kinetic energy E{x) is , see 



equation (5.69) in Pope (20001 , 



E{x) 



irp 



o^(x) 2 Jo 



(17) 



which has units of 



The previous formulas are exactly the same as in 



Pope ( 2000 1 ; we now continue toward the astrophysical 



applications. The self-similar solution for the velocity , 
equation ^ , can be re-expressed introducing the half 
width X ^ hi / S 



1 



(18) 



where A = y/2 — 1 . From the previous formula is clear 
the universal scaling of the profile in velocity that is 
reported in Figure [l] 

From a careful inspection of the previous formula 
it is clear that the variable x should be expressed in 
d units (the nozzle's diameter) in order to reproduce 
the laboratory results. In doing so we should find the 
constant k that allows us to deduce B 



B = k X SaxVT 



(19) 



Table [T] reports a set of S , B and i't for different 
opening angles a . The assumption here used is that k 
is the same for different angles. The velocity expressed 
in these practical units is 



<U >^ BUi 



x{l + A{f^r) 



(20) 



This formula can be used for x expressed in d-units 
when X > B. 

The first derivative of the profile in velocity as given 



by formula ( 20 1 with respect to the radius is 

d 



dr 



<U>-- 



Ui X 



-4 {V2-1)A kb^'^A rd 
tan(f) - I) X {b^^ + A r'^^^ bj_' 



(21) 
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The production of turbulent kinetic energy is 




Fig. 1 Mean velocity profile vs. channel radius in a tur- 
bulent round jet. The velocity distribution is a function 
growing from zero (at the wall of the channel) to a maxi- 
mum value in the central region. The experimental data at 
Re 10^ and x/d= 40 , 60, 75 and 97.5 can be found in 
Wygnanski Fiedler 1969. 



Table 1 Parameters of the turbulent jet when k — 0.54 



a[rad] 


a[degree] 


S 


B 


vt 


0.035 


2 


0.017 


30.96 


0.0052 


0.087 


5 


0.043 


12.37 


0.013 


0.14 


8 


0.07 


7.72 


0.021 


0.185 


10.64 


0.093 


5.79 


0.028 


0.261 


15 


0.131 


4.1 


0.039 


0.343 


20 


0.17 


3.06 


0.053 


0.436 


25 


0.22 


2.43 


0.067 


0.523 


30 


0.26 


2.01 


0.08 



-p = lyj^U^ X 



2 (V2-1) A k^bi'^A^^d^ 



(22) 



tan 



It is interesting to note that the maximum of T-" , is at 
1 



0.69&1 



(23) 



^/b^/A 

In these practical units the flow rate of mass , equa- 



tion (14), becomes 

TT px (tan (f ))'^ Ui d 



rh{x) = 



and flow of kinetic energy , equation (17 1 , is 



d3[/j3(tan(f))%^ 



10 



(24) 



(25) 



A more sophisticated approach makes extensive use of 
a statistical mass-averaging technique for compressible 



turbulent flow, see Favre, A. (19691; Bicknell (19841. 



3.2 The simple solution 

We now outline the conservation of the momentum flux 



in a "turbulent jet" , see Landau (19871 . The initial 



point is characterized by the following section 

Ao^nrl . (26) 

On introducing a ,the opening angle , xq ,the initial 
position on the x-axis, and vq , the initial velocity , the 
section A at position x is 



A{x) = 7r(ro -I- (x- a;o)tan (^))' 



(27) 



The conservation of the total momentum flux states 
that 



pv^Ao = pv{x)'^A{x) , 



(28) 



where v{x) is the velocity at position x . Due to the 
turbulent transfer, the density p is the same on both 



the two sides of equation ( 28 1 . The trajectory of the 



jet as a function of the time is easily deduced from 



equation (28 1 



xo tan (f ) -ro + ^ rg {rg + 2 tan (f ) vg t) 
tan(f) 



(29) 
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v{t) 



The velocity as function of the time turns out to be 

■"0 ro 



^ro (ro + 2 tan (f)?;ot) 



(30) 



The flow rate of mass and kinetic energy are respec- 
tively 



m{x) = nr'^pv , 



(31) 



(32) 



where r and v are the momentary radius and velocity 
of the jet. 



4 The physics of HH's 

This Section reports the centerline velocity, the equa- 
tion of motion , the flow of mass and the flow of energy 
for the two turbulent models here considered. 

4.1 The exact solution 

Equation (|8| allows us to deduce the centerline velocity 
of the turbulent astrophysical jet 



uo{x) 



0.5409 Mjoodj Km 
tan (0.00872 a^eg) 2;j sec 



(33) 



where adeg is the opening angle expressed in degree, 
Uioo is the initial velocity expressed in units of 100^ , 

uioo = ^^^100^^^ I '^1 is the diameter of the nozzle in pc 
units and a;i is the length of the jet in pc units. 

The previous equation allows us to deduce the equa- 
tion of motion for a turbulent astrophysical jet , 



x{t) = 1.050, 



pc 



(34) 



tan (0.00872 adeg) 
where t4 = t[yr]/{10*) . The radius of the turbulent jet 



di 
^ 2 



(35) 



Combining equations p3| , (34 1 and (pM| is possible 



to deduce the velocity of the HH object , for example 
HH34 , as function of the time see Figure [2] 
The power released in the turbulent cascade is 



£{r,x) 



2 {V2-l) kH.^^A 



12 A 2^2 



tan 



(f)(V2-l)^(6.^ + Ar2) 



x(-) 

Xi 



The flow rate of mass , see equation ( 24 ) , as expressed 
in these astrophysical units is 



rji(x) — 

0.1910 p, no xi (tan (0.008725 adeg))^ uioodi 

year4 ' 



(37) 



where no is the number density expressed in parti- 
cles cm~^ (density p ~ n-om, where m = pmYi) , p 



is the mean molecular weight (see McCray and Layzer 
(19871 suggests ^=1.4 ) , mn is the hydrogen mass , 
M.Q is the mass of the sun and yr4 are 10^ year. On in- 
troducing the solar system abundances , N{El) , where 
El represents the considered element , see Table 2 in 
and the time expressed in year we 



Lodders (2003) 



obtain 



'm{x) = 

1.9 X 10"^ p no xi (tan (0.008725 adeg)f uwo di 

N{El) Me 



N{H) year 



(38) 



where N{H) is the Hydrogen solar system abundance. 
As an example when the Fe is considered we obtain 



6.59 X 10" 



-10 



rh{x) = 

pnoxi (tan (0.008725 ttdeg))^ iijoo 



year 
case of Fe 



, (39) 



where §^ 



8.380x10° 
2.431x1010 ' 



see Table 2 in 



(20031. 



Lodders 



A comparison of the previous formula can be done 
with M{[FeII]) that in HHl varies between 2.2 x 
IQ-'^MQyr-^ and 2.8 x IQ-^^Mayr'^ , see Table 3 in 



Nisini et al. (20051 




' Time (fo**4 yr)' 



Fig. 2 Velocity of HH34 in 



versus time in t4 units 



(36) wlien Q = 2.86° , uioo = 6.5 , di = 1.94 10"^ and xi = 1.49. 
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The flow of energy , equation (251, in these astrophys- 
ical units is 

1.203 X 10^4^ no {tan {0.00872 adeg)f uioo^di^ 



ergs 



(40) 



The analysis makes extensive use of Favre's (1969) sta- 
tistical mass-averaging technique for compressible tur- 
bulent flow 

4.2 The simple solution 



Equation (|30| allows us to deduce the centerline veloc- 

ito(x) = 



ity in the simple case 

1.021uioodi 



di +2 tan (0 .00872a deg) xi - di tan (0.00872ade9) 

Km 



sec 



(41) 



The astrophysical version of the equation of mo- 

di 



tion, formula (29 1, is 



x{t) = 0.5di - 0.5 



tan (0.00872Q;deg) 



0.707 



^Jo.bdi'^ + 2.042di tan {0. 00872a deg) uiooU 
tan(0.00872a<ie3) 

The flow rate of mass , see equation (pTl) , is 



(42) 



m{x) = 0.01259 (0.5 di + tan (0.00872 adeg) xi 



-0.5 dj tan (0.00872 a(ie3))(^j uioo i^f^o 



Me 
year4 



The flow of kinetic energy ,equation (|32| , is 

E{x)^ 



6.231 X 10^3^710 uioo^di^ 



ergs 



di + 2.0 tan (0.00872 Udeg) xi 



(43) 



(44) 



5 Complex trajectories 

This section reports the kinematic effects that lead to 
complicate trajectories as well an explanation for the 
train of knots which are visible in the first part of the 
HH objects. 

5.1 The precessing jets 

The wide spectrum of observed morphologies that char- 
acterizes the HH objects can be due to the kinematic 



effects as given by the composition of the velocities of 
different effects such as decreasing jet velocity , jet pre- 
cession and proper velocity of the host star in the in- 
terstellar medium (ISM). Of particular interest is the 
evaluation of various matrices that will enable us to 
cause a transformation from the inertial coordinate sys- 
tem of the jet to the coordinate system in which the 
host star is moving in space. The various coordinate 



systems will be x=(a;,j/, z) , ic^ '={x^ \y 



r-(l) 7,(1) ^(1) 



x(^)=(a;('^', y^'^), z^'^)). The vector representing the mo- 
tion of the jet is represented by the following 1x3 
matrix 



x{t) 



(45) 



where the jet motion L(t) is considered along the x-axis. 

The jet axis, x, is inclined at an angle ^prec relative 
to an axis x^^^ and therefore the 3x3 matrix, repre- 
senting a rotation through the z axis, is given by: 



F = 



cos {^p 

rec 

) -Sm {^prec) 
sin {^prec) COS prec) 











(46) 



From a practical point of view ^prec can be derived by 
measuring the half opening angle of the maximum of 
the sinusoidal oscillations that characterizes the jet. 

If the jet is undergoing precession around the x^^^ 
axis, riprec can be the angular velocity of precession 
expressed in radians per unit time ; ilp^ec is computed 
from the optical maps by measuring the number of sinu- 
soidal oscillations that characterize the jet. The trans- 
formation from the coordinates x'^) fixed in the frame 
of the precessing jet to the non-precessing coordinate 
x'-^) is represented by the 3x3 matrix 



P 











COS iyiprect) — siu (fiprec^) 
sin {^prect) COS (f2prec^) 



(47) 



As an example Figure |3] reports the precessing jet ap- 
phed to HH34. 

The last translation represents the change of frame- 
work from (x^'^-'), which is co-moving with the host star, 
to a system (x^^)) in comparison to which the host star 
is in a uniform motion. In the laboratory experiments 
the velocity of the host star is replaced by a wind , see 
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Figure 3 in Ciardi et al. (20081. The relative motion 
of the origin of the coordinate system {x^^\ y*^'^^ z*^^^) is 
defined by the Cartesian components of the star veloc- 
ity Vx,Vy,Vz, and the required 1x3 matrix transforma- 
tion representing this translation is: 



B 



Vrr.t 



Vyt 



V,t 



(48) 



On assuming, for the sake of simplicity, that Vx~0 and 
Vz^O, the translation matrix becomes: 



B = 







Vyt 







(49) 



In other words, the direction of the star motion in the 
ISM and the direction of the jet are perpendicular. 
From a practical point of view the star velocity can 
be measured by dividing the length of the star in a di- 
rection perpendicular to the initial jet velocity by the 
lifetime of the jet. The final 1x3 matrix A representing 
the "motion law" can be found by composing the four 
matrices already described 

A = B + {P-F-G) 



cos i^prec) X (t) 
Vyt + COS (ilprect) siu prec) X (t) 
sin (flprect) sin {^prec) X (t) 



(50) 



The three components of the previous 1x3 matrix A 
represent the jet motion along the Cartesian coordi- 
nates as given by the observer that sees the star moving 
in a uniform motion. As an example Figure |4] reports 
the effect of inserting the star's velocity on the precess- 
ing IIH34 as plotted in Figure |3] 

The fifth matrix allows to model the point of view 
of the observer through the matrix E representing the 
three Eulerian angles which characterizes the point of 
view of the observer, 6,<i>,^' , see [Goldstein, Poole,] 
and Safko| ( |2002[ ) . The product E ■ A is not reported 



for space problem and Figure |5] reports the same as 
Figure |4], but from a particular point of view. In other 
words the particular point of view can produce complex 
projected patterns of a simple basic trajectory as rep- 
resented by Figure |4j A comparison of Figure |5] should 
be done with the image of HH34 as available at 
|http: / / antwrp.gsfc.nasa.gov/ apod/ap991129.html| 




Fig. 3 Continuous three-dimensional trajectory of HH34: 
the three Eulerian angles characterizing the point of view 
are $= ° , Q— ° and °. The precession is charac- 
terized by the angle 'i'prec= 10 ° and by the angular velocity 
36.00 [°/10*year]. The physical parameters charac- 



terizing the jet motion are : itioo 
di = 9.7 10-* and qo= 2.86 °. 



: 6.5, t4=4 ,xi = 1.49 



minrmum= 0.0556 maximum= 1.0000 




[pc] 



Fig. 4 Continuous three-dimensional trajectory of HH34: 
the three Eulerian angles characterizing the point of view 
are ° , 6= ° and °. The precession is 
characterized by the angle ^'prec= 10 ° and by the angu- 
lar velocity Qprec— 36.00 [°/10*year]. The star has veloc- 
ity Vy — 11.19^^^. The physical parameters characteriz- 
ing the jet motion are : uioo ~ 6.5, ti = 4 , xi = 1.49 , 
di = 9.7 10"* and ao= 2.86 °. Here the plane of the trajec- 
tory is perpendicular to the observer . The two Cartesian 
axis are in pc and the integral operation which allows to 
build the image is performed on cubic grid of 1200^ pixels. 
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made with the VLT by the FORS Team or Figure 1 
in Reipurth et al. (2002) which has a field of 4.5' x 5'. 



The astrophysical version of the star's motion as rep- 



resented by the translation matrix B, formula (49 1, is 



These observable quantities are identified with the 
following theoretical variables: 



(53) 



y = 0.01021 Vy t4 , (51) 
where Vy is expressed in units and t4 = t[yr] / (W^) 

The previous equation can be combined with the mo- 



tion along X as represented by equation (34) in order 
to find the angle /3 in degree that characterizes the tra- 
jectory: 



57.29 arctan( 



13- 

0.009718^;y^/^ 



arctan( — ) 

X 



y/UiooVd^^/tB,n (0.008727adeg) 



) °- (52) 



This angle varies from when t4—0 to 23.39° when t4—4 
and the parameters of Figure ^ are used. 

Is also interesting to point out that a rotation of 
90° around the y axis of the trajectory as reported in 
Figure |4] makes the jet straight rather than bended. 

Analogous results on ballistic jets from precessing 
sources has been obtained, see |Lightfoot and Glencross] 
([T986|;lRaia|p993|. 



5.2 The Kelvin-Helmholtz instabilities 

The macroscopic phenomena of the jets as the pres- 
ence of knots and wiggles can be due to the Kelvin- 



Helmholtz instability (after Kelvin (1871); Helmholtz 



(1868 1 ) of an axisymmetric flow along the velocity-axis 
when the wavelengths A = ^ ( fc is the wave-vector) 
are greater than the jet radius Uj, which is taken to 



be independent of the position along the jet, see Fer- 



rari, Trussoni, and Zaninetti 


(1979 


1981 


I; 


Ray and Er- 


shkovich 


( 


1983 


); 


Hardee, Clarke, and Howell 


( 


1995 


). 



The velocity , Uq, is assumed to be rectangular. The 
internal ( external ) fluid density is represented by poi 
(poe) , the internal sound velocity is Si and vo= 
Starting from the equations of motion and continuity, 
and assuming both fluids to be adiabatically compress- 
ible, it is possible to derive and to solve the dispersion 
relation from a numerical point of view , see |Zaninetti| 



(19871 



We then start from observable quantities that can be 
measured on radio-maps such as the total length Lo^s, 
the wavelength A^;,^ of the wiggles (m=l) along the jet, 
the distance A^j,^ (m=0) between knots, and the final 
offset A Lobs of the center of the jet. 



^7nax ^obs 



^0 exp{— — -) 



'MtadCij ' 



AL 



obs 



flip — Lf-, 



(54) 



(55) 



(56) 



where tad — imin • Si I Uj and Aq is the amplitude of the 
perturbed energy. The result is a theoretical expres- 
sion for tmin the minimum time scale of the instability, 



the wavelength connected with the most unstable 



An 

mode and le the distance over which the most unstable 

These 



mode grows by a factor e, see Zaninetti (1987). 
parameters can then be found through the set of non- 
linear equations previously reported. By choosing two 
objects, HHl and HH34 the observational parameters 
can be measured on the optical image, see Table |2] 

The four nonlinear equations are then solved and the 
four theoretical parameters are found , see Table |3] 

An application of the results for HHl here obtained 
is reported in Figure [6] ; the comparison should be done 
{[SI I]) in [Reipurth erd] (|2000[) that 



1 



with Figure 

covers « 14.16 arcseconds. The application to HH34 is 
reported in Figur e [7] and the compa rison should be done 
with Figure 3 in Reipurth et al. (2002) which covers 



« 30 arcsec. In both cases the wavelength of the pinch 
modes (m = 0) and the oscillations of the helical mode 
(m = l)are those reported in Table [3] 



Table 2 Parameters of the observed oscillations in HH 
objects expressed in averaged radius units 



Geometrical measure 


HHl 


HHM 


Lobs[o,veraged radius units] 


73.3 


54 


Aobs [averaged radius units] 


8.1 


4.5 


Xlij^[averaged radius units] 


29.3 


27 


ALobs[o-veraged radius units] 


1 


1 
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minrmum= 



0.0098 



maxirrium = 



1.0000 




I 




'x [radius units J" 



Fig. 5 Continuous three-dimensional trajectory of HH34: 
the three Eulerian angles characterizing the point of view 
are <E>= 290 ° , 9= 180 ° and 18 °. The precession 
is characterized by the angle 'i/prec— 10 ° and by the an- 
gular velocity flprec= 36.00 ["/IC^year]. The star has ve- 
locity Vy — 11.19^^^. The physical parameters character- 
izing the jet motion are : uioo = 6.5, t4 = 4 , xi = 1.49 , 
di = 9.7 10-* and qo= 2.86 °. The image is here reported 
as a theoretical 2D surface brightness of emission, the two 
Cartesian axis are in pc and the integral operation is per- 
formed on cubic grid of 1200^ pixels. 



Fig. 6 Superposition of the pinch mode (m = 0) on the 
helical mode (m — 1) for HHl in radius units. The param- 
eters are as in Tables [2] and [S] 



minimurri= 0.001 ' 



ma>;imum= 0.2318 



Table 3 Theoretical parameters from oscillations deduced 
from the four nonlinear equations 



Theoretical variable HHl HHiA 



n 
M 



1.22 
52.2 

325 



1.96 
26.67 
15.67 




I 



20 25 



0.012 0.017 



Fig. 7 Superposition of the pinch mode (m = 0) on the 
helical mode (m = 1) for HH34 in radius units. The param- 
eters are as in Tables [2] and |3] The image is here reported as 
a theoretical 2D surface brightness of emission, the integral 
operation is performed on cubic grid of 1200^ pixels. The 
intensity or brightness in the spherical blobs (m = 0) scales 
as in equation (71 1. 



10 



6 The image from geometry 

The transfer equation in the presence of emission only 
, see for example equation (1.27) in Rybicki and Light 



man 



(2003 



1991 ) or equation (4.9) in Dopita and Sutherland 



IS 



ds 



(57) 



where Ii, is the specific intensity or brightness which 

has units of ^ ^^^J^r Hz ' ^ '^^ li'^^ '^ig^* ' ^'^ 
the emission coefficient which has units of — tj- 

s cm^ ster Hz 

, ku a mass absorption coefficient, C the mass density 
at position s and the index v denotes the interested 



frequency of emission. The solution to equation (57 1 is 



/.(r.) = f (l-e--'^(^)) , 

where Ti, is the optical depth at frequency u 
dTi, — k^Cds . 



(58) 



(59) 



We now continue analyzing the case of an optically thin 
layer in which Ti, is very small ( or k^, very small ) and 
the density C is substituted with our number density 
C(s) of particles. Two cases are taken into account : 
the emissivity is proportional to the number density 
and the emissivity is proportional to the square of the 
number density . In the linear case 



where K is a, constant function. 
In the quadratic case 



(60) 



where K2 is a constant function. This is true for ex- 
ample for free-free radiation from a thermal plasma, 



see formula (1.219) in Lang (19991 or formula (6.17) in 
Dopita and Sutherland ( 2003 1 . 



The intensity is now 

I^{s) I C{si)dsi 

J So 

optically thin layer linear case , 



(62) 



or 



his) = K2 f C{sifdsi 

J So 

optically thin layer quadratic case 



(63) 



In the Monte Carlo experiments the number density is 
memorized on a 3D grid M(i,j, k) where i,j and k are 
indexes varying from 1 to pixels , and the intensity is 



I{i,j)=Y,tis-KM{i,],k) 

k 

optically thin layer linear case , 



(64) 



or 



/(z,j) = ^AsxA^(j,j,fc)2 

k 

optically thin layer quadratic case , 



(65) 



where As is the spatial interval between the various 
values of intensity and the sum is performed over the 
interval of existence of the index k. In this grid frame- 
work the little squares that characterized by the posi- 
tion of the indexes i,j correspond to a different line of 
sight. When all the different pixels are viewed together 
the image is formed. The ensemble of all the pixels can 
be considered a theoretical surface intensity or a the- 
oretical surface brightness. We now outline a possible 
source of radiation. The volume emission coefficient of 
the transition j2i is 



J21 



712^21^1^21 

47r 



(66) 



where level 1 is the lower level , level 2 is the upper level 
, 712 is gas number density , 712^21 the rate of photons 
emitted from a unit volume , A21 is the Einstein coef- 
ficient for the transition , h is the Planck constant and 



V21 the considered frequency, see Hartigan (20081. In 



the case of optically thin medium the intensity of the 
emission I21 is the integral along the line of sight 



(61) /21 - / hidl 



In the case of constant gas number density 
/21 oc Z , 



(67) 



(68) 



where I is the considered length that in the astrophysi- 
cal diffuse objects depends from the point of view of the 
observer. The optically thin layer approximation repre- 
sents therefore a useful approximation to build models 
for the intensity of radiation which are comparable to 
the observed profiles. 

We now analyze the behavior of the intensity of a 
cross section of a jet , the behavior of the maximum 
intensity at the centerline as a function of the distance 
from the central source , the intensity of complex mor- 
phologies and the sudden increase in intensity as given 
by the toroidal jet. 
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6.1 Intensity at a fixed distance 

We explore the behavior of the intensity or brightness 
along the a jet when the distance from the origin ,y, 
is fixed . We assume that the number density C is 
constant in a cross section of radius a and then falls to 
, see Figure |8] 

The length of sight , when the observer is situated 
at the infinity of the x-axis , is the locus parallel to the 
X-axis which crosses the position ?/ in a Cartesian x — y 
plane and terminates at the external circle of radius a. 
The locus length is 



(69) 



When the number density C,„ is constant in the cylin- 
der of radius a the intensity or brightness of radiation 
is 



/oa = C™ X 2 X (Va'-z/2) ; < 2/ < 



ha = Cm X 2 X a X cos(a) ; -| < a < | 



(70) 



(71) 



that can be named the "trigonometrical law" for the 
intensity or brightness . Is interesting to underline that 
the two previous equations hold for a cylindrical and a 
conical jet as well for a spherical blob when the number 
density is constant. 

6.2 Centerline Intensity function of the distance 

We now explore the behavior of the intensity or bright- 
ness at the centerline of the jet as a function of the 
distance x from the nozzle. From the previous para- 
graph |6.1| we learned that the maximum intensity or 
brightness at the centerline of the jet at a fixed dis- 
tance X is proportional , as a first parameter , to the 
jet's diameter d , 



I{r = 0, a;) cx d{x)Cm oc xC„ 



(72) 



In order to have a constant intensity or brightness along 
the centerline of the jet as function of x , the number 
density Cm of the emitting particles should decreases 
as 



1 



Cm{x) cx - 

X 

As a consequence the intensity or brightness 



l[x) oc Lm y. d (X X , 

X 



(73) 



(74) 



will be constant along the jet. In the framework of 
the optically thin medium the emitting length will not 



change but conversely the number density can take the 
general form 



1 



Cmix) CX 



(75) 



which means that the intensity or brightness scales as 



1 



I{x) cx d{x)—^^ cx X 



(76) 



The value of s can found from the scaling of the ob- 
served intensity or brightness as function of x. As an 
example of constant intensity or brightness of emission 
along a knotty jet we report the image of the first part 
of HH34 where Cm{x) cx 1^0.15 was chosen ; a com- 



parison should be done with Figure 3 in |Reipurth et al\ 
(20021. The relative cut along the jet's main axis of 
symmetry, is reported in Figure [9] where s = 0.15 was 
used. As a practical example the intensity (counts) of 
He for HHllO as in Figure 3 of |Riera et aI1 ( p003| de- 
creases of a factor w 8.6 from the second blob B to the 
last blob Q. In our simulation of HH34 as reported in 
Figure |9] the intensity decreases of a factor « 8.3 from 
the first to the last blob. 

The 3D algorithm already presented replaces previ- 
ous efforts based on the 2D random walk , see Figure 



7 in Zaninetti ( 1989 1 and on the 3D random walk from 



many injection points , see Figure 8 in |Zaninetti| ( |1999| ). 
6.3 Complex Morphologies 

The integral operation of the emissivity along the line 
of sight of a turbulent jet can be performed in an an- 
alytical way only in a simple configuration : the jet 
perpendicular to the observer, see Section [7?T| The con- 
currency of complex trajectories and a general point of 
view of the observer characterized by the three Eulerian 
angles Q, <i> and asks a numerical treatment. We re- 
member that the points that characterize the trajectory 
of HH34 , see Section |5.1| , are already in such a way 
that the product e(x) xx is nearly constant. This means 
that the intensity or brightness is nearly constant along 
the main direction. These points are inserted on a 3D 
grid made by pixels^ points and a sum is performed 
over one index, see Figure [5] 

The enhancement of the intensity or brightness in 
the previous map where the jet is bending is due to 



the particular point of view of the observer. Figure 10 
reports a cut along the centerline of a jet from which 
is possible to observe an increase of a factor « 5 in 
the axial intensity or brightness otherwise constant . 
An analytical evaluation of such increase is reported in 
Section iMl 
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Fig. 8 The source is represented through a circular section 
perpendicular to the jet axis. The observer is situated along 
the X direction, one line of sight is indicated and the angle 
a is clearly indicated. 




X tarcsec] 



Fig. 9 Intensity of HH34 represented through a cut in the 
x-direction (y=0) , parameters as in Figure]?] The intensity 
or brightness along the cut decreases as I{x) oc x~°'^^. 




x°fpc] 



Fig. 10 Behavior of the intensity or brightness of HH34 
at the centerline. Parameters as in Figure [5] 



Conversely when the plane of the trajectory is per- 
pendicular to the point of view of the observer the en- 
hancement in the intensity or brightness of HH34 is not 
present , see Figure |4] 

6.4 Toroidal Model 

The curved shape of a jet of finite cross section is not 
easy to parametrize. The torus represents a possible 
model due to the presence of the small radius that char- 
acterizes the cross section of the HH object , r , and the 
great radius R that can be identified with the curvature 
k, k = 1/R that characterizes the 3D trajectory. The 
torus has the following parametric equations: 



X — cos(s) ■ {R+r ■ cos(i)) 
y — sin(s) ■ {R + r ■ cos(t)) 
z — r ■ sin(t) , 



(77) 



where t g [0,27r) and s G [0, 7r/4). 

Figure [TT] reports a section in the middle of the torus 
z = 0, from which is possible to see that the dotted 
line presents the longest line of sight , Imax , when the 
observer is at infinity of the x — axis . The shortest 
line of sight is 2r . The maximum enhancement in the 
presence of constant number density , e , is 



2r 



A simple geometrical demonstration gives 



(78) 



(79) 



The radius that produces an enhancement e in the in- 
tensity or brightness is therefore 



R 

r 



(80) 



As an example an enhancement of e = 5 is produced 
by a radius of curvature 25 times greater in respect to 
the HH's radius. 



7 The image from turbulence 

The power released in the turbulent cascade has the 
same dimension of the emission coefficient and therefore 
can be considered the source of emissivity. We now 
consider a linear and a nonlinear relationship between 
turbulent power and emission coefficient. 
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7.1 Linear correspondence 



It is assumed that the emission coefficient of the HH 
scales as the power released in turbulent kinetic energy, 



see equation ( 22 1 



(81) 



Due to the additive property of the optically thin 
medium along the line of sight, an integral operation is 
performed in order to obtain the intensity or brightness 
of emission 



I{y) ^ I 2 X e{r)dz 

Jo 



(82) 



with r = -^/z^ + and a representing the jet radius, 
see Figure |8] 

The intensity or brightness of emission according to 
formula (1221) is 



liy) 



Il{z)dz (83) 



with Il{z) 



(hi.^ + Az'^ + Ay-^'^ x2 (\/2- l)tan(f) 

This integral has an analytical solution but it is com- 
plicated and therefore Figure [12] only shows the numeri- 
cal integration which presents a characteristic shape on 
the top of the blob called the " valley on the top" . The 
maximum of this integral is at the point y « 0.495i 
and the value of intensity or brightness at the maxi- 
mum is 1.18 times the value at the point ?/ = ( the 
center of the jet). The near infrared images of HH 110 
jet were interpreted as due to low velocity shocks pro- 
duced by turbulent processes, see |Noriega-Crespo et al.\ 
(19961. The spatial intensity or brightness distribution 
of H2 , Ha and [^z/Jgyiy/si perpendicular to the flow 
axis and along the cross section of knots in HHWi) has 
a behavior that can be approximated by a Gaussian dis- 
tribution , see Figure 4 in ,Noriega-Crespo et al.^ (1996[ ) . 
In one case , Ha in knot P of HHllO in Figure 5 in 
Noriega-Crespo et al. (19961 , it is possible to see a 



bump near the maximum of the intensity or brightness 
in the transversal direction. A second observation that 
presents a bimodal profile is the spatial Ha intensity 
distribution through the cross section of knot / -I- J of 
HHllO visible in Figure 10 in [Riera et al\ ( |2003[ ). A 
third observations is the \SII] + [Nil] profile in blob 1 
of HHllO as in Figure 17 oflHartigan et aZ.|(|2009l. 



These three cases can be considered an observational 
evidence of the physical effect previously named " valley 
on the top" . 




Fig. 11 The section of one fourth of a torus is represented 
through a circle of radius R and a bigger circle of radius 
R -\- 2r. The observer is situated along the x direction, the 
line of sight of maximum length is indicated. 




y [b,/2] 



Fig. 12 Intensity of radiation I{y) in 61 units across the 
turbulent jet when z/d = 50 , k =0.54, A = 0.414 and 
ctdeg = 2.0 : full line ( linear correspondence) and dashed 
line (non linear correspondence,/ — 4) . 
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Is also possible to build a 2D map of the surface 
brightness of emission computed according to the inte- 



gral of equation! 83) in a conical jet as HHl and Fig- 



ure [13] reports such a map. 

In this case the values of emissivity are memorized 
on a 3D grid made by (pixels)^ points. The integral is 
represented by a sum along the line of sight and Fig- 
ure [14] reports a cut in the middle. 

7.2 Non linear correspondence 



From a careful inspection of formula ( 22 ) we see that 



the local power released in the turbulent cascade at 
r — bi scale as ~ . In order to have a constant 
intensity or brightness along the jet we now consider 
the case 

e^iVy/^ , (84) 
which means that the intensity or brightness scales as 



X 

The intensity or brightness of emission is 

Jo 



(85) 



(4 



In{z)dz (86) 
with In{z) = 

)^. 



22 + A 2;2 (V2- l)tan(f) 



The numerical result for f ~ A , constant intensity or 
brightness along the main direction x, is reported in 



Figure 12 as a dashed line. 



8 Conclusions 

Law of motion The two theories here considered are 
based on the behavior of the centerline velocity , see 
the astrophysical equations (33 1 and (41 1 . From the 



two previous equations is possible to deduce the law of 
motion in presence of a stationary state, see the astro- 
physical equations ( 34 1 and ( 42 1 . On the way the flow 



rate of mass and the flow rate of energy (the mechani- 
cal luminosity) are also derived , see the astrophysical 
equations (pfl ,(l43| , and (liil. 



Images The analysis of the intensity or brightness 
of a HH object has been split in three theoretical parts 
corresponding to three observable cases. 

1. Transversal cut We have analyzed the case of con- 
stant number density , see equation (71), and the 



minimum^ 



0.0047 



ma>!imum = 



1.0000 




I 
I 



0.5 1 1.5 

z [radius units] 



Fig. 13 Theoretical 2D map of the surface brightness 
of emission corresponding to a 2.2 arcsec of HHl around 
z/d ^ 110. The parameters are k =0.54, A = 0.414 and 
ctdeg ~ 2.08 . The integral operation is performed on a 
cubic grid of 400^ pixels. 




arcsec 



Fig. 14 Intensity of HHl represented through a cut in the 
y-direction , parameters as in Figure [l3] 
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case of emissivity connected with the power released 
in the turbulent cascade, see equations ( 83 ) and ( 86 ). 
The intensity or brightness from turbulent cascade 
originates a curious effect at the center of the jet 
named "valley on the top" . 

Longitudinal cut Through a parametrization of the 
number density is possible to fit the theoretical and 
the observed intensity or brightness , see equation 
(|76| . 

2D map The details of the HH's image can be sim- 
ulated imposing an arbitrary point of view of the 
observer. The enhancement in intensity or bright- 
ness is explored from a numerical point of view , 
see Figure ^ and (10). An analytical explanation 
of the enhancement in intensity or brightness is de- 
rived from the geometrical properties of the torus 
, see formula (79). Is interesting to underline that 
the " torus effect" replaces the concept of bow shock, 
that is often used in order to explain the intensity 
enhancement along the HHs, see Krist et al. (1999); 



Smith, Khanzadyan, and Davis (20031. 
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